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The question of whether one should speak of a "pairing glue" in the Hubbard and t-J models is 
basically a question about the dynamics of the pairing interaction. If the dynamics of the pairing 
interaction arises from virtual states, whose energies correspond to the Mott gap, and give rise to 
the exchange coupling J, the interaction is instantaneous on the relative time scales of interest. 
In this case, while one might speak of an "instantaneous glue", this interaction differs from the 
traditional picture of a retarded pairing interaction. However, if the energies correspond to the 
spectrum seen in the dynamic spin susceptibility, then the interaction is retarded and one speaks of 
a spin-fluctuation glue which mediates the d-wave pairing. Here we present results from numerical 
studies which provide insight into this question. 



The question of whether the pairing interaction in the 
cuprate superconductors should be characterized as aris- 
ing from a "pairing glue" has recently been raised (ij. 
As we will discuss, this is a question about the dy- 
namics of the pairing interaction and it will be an- 
swered when we know more about the frequency depen- 
dence of the cuprate superconducting gap. From the 
d-wave (cos kx cos ky) momentum dependence of the 
cuprate gap, we know that the pairing interaction is spa- 
tially a short range, dominantly near- neighbor attraction. 
However, in spite of pioneering ARPES[iJ|,Jl i, 0, 
tunneling 0, and infrared conductivity [l3, EH stud- 
ies, we do not yet have sufficient information to defini- 
tively characterize its dynamics. Thus, while there is 
a growing consensus that superconductivity in the high 
Tc cuprates arises from strong short-range Coulomb in- 
teractions between electrons rather than the traditional 
electron-phonon interaction, the precise nature of the 
pairing interaction remains controversial. 

This is the case even among those who agree that the 
essential physics of the cuprates is contained in the Hub- 
bard and t-J models. For example, both Anderson's 
resonating-valence-bond (RVB) the orvll2l| and the spin- 
fluctuation exchange theory[lJ, [3, Il5l| lead to a short- 
range interaction which forms dx^-y^ pairs. However, 
the dynamics of the two interactions differ. In the RVB 
picture, the superconducting phase is envisioned as aris- 
ing out of a Mott-liquid of singlet pairs. These pairs are 
bound by a superexchange interaction J which is propor- 
tional to t'^ /U. Here t is the effective hopping matrix ele- 
ment between adjacent sites and U is an onsite Coulomb 
interaction. J is determined by the virtual hopping of an 
electron of a given spin to an adjacent site containing an 
electron with an opposite spin. [16| Thus the dynamics of 
J involves virtual excitations above the Mott gap which 
is set by U, and the pairing interaction is essentially in- 



stantaneous. In this Anderson has noted, [l| one 

should not speak of a "pairing glue" in the same sense 
that this term is used when referring to a phonon me- 
diated interaction. In the spin-fluctuation exchange pic- 
ture, the pairing is viewed as arising from the exchange of 
particle-hole spin fluctuations whose dynamics reflect the 
frequency spectrum seen in inelastic magnetic neutron 
scattering. This spectrum covers an energy range which 
is small compared with U or the bare bandwidth 8t. In 
this case, the pairing interaction is retarded and in anal- 
ogy to the traditional phonon mediated pairing, one says 
that the spin- fluctuations provide the "pairing glue." So 
the question of whether there is a "pairing glue" offers a 
way of distinguishing different pairing mechanisms. Here 
using numerical techniques we examine this question for 
the t-J and Hubbard models. 

In the superconducting state the Nambu self-energy 
Y,{k,LL!) can be parametrized as 

E(A:,a;) = (1 - Z{k,uj))ujTQ -\- x{k,uo)Ts, + (f>{k,uj)Ti. (1) 

Here, tq, ti and T3 are the Pauli spin matrices and Z{k, lu) 
and xikjUj) describe the so-called normal components of 
the self-energy and the gap function (l){k,uj) describes 
the anomalous part which contains information on the 
internal structure of the pairs. The complex gap func- 
tion (f>(k,u!) = 4'i{k,Lu) -\- i(j)2{k,Lu) satisfies the Cauchy 
relation 

Mk.^)^'-l^-^dJ (2) 
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and for = 0, one has 
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Based upon this, a useful measure of the frequency de- 
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pendence of the pairing interaction 17| is 

2 02(fc,cj') ^_/ 
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It gives the fraction of the zero frequency gap function 
which arises from frequencies below f2. In order to obtain 
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FIG. 1: (a) The imaginary part of the Pb gap function <i>i(u)) 
versus tj (solid curve). The peaks in t^-iiyi) occur at the trans- 
verse and longitudinal ojl peaks of a^_F(tj) (dashed curve) 
shifted up by the gap Aq. (b) The pairing interaction spec- 
tral weight /(f2) versus for Pb. liVC) increases as passes 
through ojt + Ao and wz, + Aq reflecting the transverse and 
longitudinal phonon contributions to the pairing. At larger 
values of fi, /(SI) exceeds unity because 01(0) is reduced from 
the value it would have just due to the phonons by the pres- 
ence of the non-retarded screened Coulomb pseudo-potential 



hopping t. The Hubbard model has an onsite Coulomb 
interaction V and its Hamiltonian is 

H -i ^(c+Cjs + c+c,^) + C/^ni|n,x-^^n,^ (6) 

(ij)s i is 

with fi a chemical potential which sets the site filling (n) . 
Here creates an electron of spin s on site i and riis — 
c^Cis is the site occupation number operator for spin s. 
The t-J model is the large U limit of the Hubbard model 
in which no double site occupancy is allowed and near 
neighbor spins are coupled by an exchange interaction J. 



H = -t 



Here S, 
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J^(S, •Sj-in.nj) (7) 



{ij) 



and is a projected fermion op- 
erator defined as c^(l — Ui-s)- 

Exact diagonalization calculations were carried out for 
the t-J model on a square cluster of = 32 sites. This 
particular cluster exhibits the full local symmetries of the 
underlying square lattice and has all of the most symmet- 
ric k points in reciprocal space. Here we will consider the 
0-, 1- and 2-hole sectors. One hole doped on to a 32-site 
cluster corresponds to a doping x ~ 0.03. 

The gap fmiction (j){k, lu) can be extracted by combin- 
ing Lanczos results for the one-electron Green's function 
G(k,uj) and Gorkov's off-diagonal Green's function [17. 



some insight into I{k, 51), we first consider the case of Pb. 
Here the k dependence of the gap function is negligible 
and only the frequency dependence enters. The imag- 
inary part of the gap function (j}2{io), determined from 
tunneling data[3| is shown as the solid curve in Fig. [T] 
The dashed curve shows a'^F{uj). Using this result for 
(/)2(w) along with the value of 0i(a; = 0), we have eval- 
uated /(51). As seen in Fig. [Hb), /(il) increases as fl 
passes through the characteristic transverse and longitu- 
dinal Pb phonon frequencies plus Aq. It then exhibits a 
broad maximum and settles down to a value that exceeds 
1. The maximum arises from the change in sign of (f>2i^) 
which occurs at a frequency 2 to 3 times the character- 
istic frequencies of the retarded part of the interaction. 
The reason that the asymptotic value of I{ft) exceeds 
unity is that the non-retarded screened Coulomb pseu- 
dopotential leads to a negative, frequency independent, 
contribution </)nr to the real part of (/)i(cj). In this case 
the Cauchy relation Eq. [5] becomes 



duj' + 0NR 



(5) 



and at high frequencies /(fl) exceeds 1 by the non- 
retarded contribution — (/'nr/</'i(0). 

The models that we will consider have a square two- 
dimensional lattice with a near neighbor one electron 



F{k,uj) = F{k,LU + iri) + F{k,LU -ir]) (8) 

with 

Fik, z) = {^oiN - 2)1 5_fc,_, -i— dka l^foiN)) 



H + En^i 

(9) 

Here the number of electrons in the initial and final 
groundstates differ by 2 and Ejss-i is defined as E^-i — 
{Eo{N) + Eo{N-2))/2. For a finite cluster, the diagonal 
Green's function is defined as 



G{k,u;) = 

{^oiN -2)\iikc 

+ {MN)\~cl 
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bj + i?] — H + En^i 

^ ~Cka\^a{N)) (10) 



uj — irj + H — E^^i 

With this definition, both G{k,uj) and F{k,uj) have the 
same set of energy poles. Using a continued fraction 
Lanczos based method both G{k,uj) and F{k,uj) have 
been calculated and the gap function (l)(k,uj) determined 
from 



(l){k,uj) = -- 



F{k,u) 



(11) 



G{k,uj)G{k,-Lo) + F^{k,uj) ■ 
Results for 0(fc, oj) and /(fc, 51) for J /t = 0.3 and x ^ 3% 
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FIG. 2; (a) The real (solid) and imaginary (dashed) parts 
of (j>2{k,Lj) versus u/t obtained for a 32-site cluster. Here 
k = (0,7r), J/t = 0.3 and the doping x ~ 3%. (b) I{k,n) 
versus Q/t for J/t = 0.3 (solid) and J/t = 0.5 (dashed) for 
k = (0, tt) and x ~ 3%. 



{N = 32) with k = (0, tt) are plotted in Fig. H We be- 
lieve that finite size effects are responsible for (f>2{k,uj) 
starting out negatively and that the corresponding neg- 
ative dip in I{k,n,) is an artifact. For J/t — 0.3, the 
rapid increase in I{k,n) as ft/t exceeds ^ 0.75 reflects 
the dynamic contributions of the spin-fluctuations and 
the broad maximum arises from the negative swing in 
(j)2{k,uj) which occurs when lu exceeds several times their 
spectral range. This is similar to the behavior seen in 
Pb when ui exceeds the spectral range of a'^F{uj). At 
higher frequencies in Fig.HlD, is seen to decrease be- 
low 1. This high frequency behavior in which drops 
below 1 is more clearly seen for J/t = 0.5, as shown by 
the dashed curve in Fig. The fact that at high fre- 
quency I{k — (0, tt), ri) lays below 1 means that there is a 
non-retarded (instantaneous) contribution to the d-wave 
pairing interaction. In contrast to the case of the tra- 
ditional low temperature superconductors, here the non- 
retarded contribution increases the pairing, correspond- 
ing to a positive value of (/)NR(A:)/0(fc, 0). Furthermore 
its relative contribution increases as J/t increases. 

For the Hubbard model, one can explore the full dy- 
namic range including the upper Hubbard band so that 
the Cauchy relation does not have an additional constant 
term (/)nr- To calculate (j>{k,uj) for the Hubbard model 



we have used a dynamic cluster approximation [20l [21 1 
(DC A). The general idea of the DCA is to approximate 
the effects of correlations in the bulk lattice with those 
on a finite size cluster with Nc sites and periodic bound- 
ary conditions. The DCA maps the bulk {L x L with 
L — > oo) lattice problem onto an effective periodic cluster 
embedded in a self-consistent dynamic mean-field that 
is designed to represent the remaining degrees of free- 
dom. The hybridization of the cluster to the host ac- 
counts for fluctuations arising from coupling between the 



cluster and the rest of the system. Here we have used 

m (NCA) to deter- 



21 



a non-crossmg approximation 
mine 4>{kA,<jj) for a 4-site 2x2 cluster at a wave vector 
kA = (0,7r). This cluster allows for a gap with d-wave 
symmetry and is such that within the non-crossing ap- 
proximation dynamic results can be obtained on the real 
frequency axis. Similar calculations were performed for 
the t-J model in Ref. 
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In mean-field theories such as the DCA, the mean-field 
generates a constant real term (pMpikA)- In the infinite 
cluster size limit, the DCA recovers the exact result and 
the mean- field contribution (pMF{kA) vanishes. For a fi- 
nite cluster size, we therefore view this contribution as an 
artifact and subtract it off of (/)i(fcAjw) before perform- 
ing the analysis based on the Cauchy relation, ^mf(^a) 
was determined from lim^^oa (f'ikAi^)- The expression 
for becomes 



I{kA,n) 



TT J{J ur TT J(J LU' 

2 TOO <hih^du^' 

TT JU CJ 



0l(fcA,O) - (pMFikA) ' 

(12) 

Results showing /(A;^,i7) versus 51 for a filling (n) = 
0.8 and U/t = 10 are plotted in Fig. [3^. The d-wave 
projection of the dynamic spin susceptibility 
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FIG. 3: (a) I{kA,^) versus Q/t for the 2 x 2 DCA-NCA 
Hubbard calculation with U/t = 10, (n) = 0.8 and T/T^ ~ 
0.95. Here kA = (0,7r). (b) The d-wave projected Xd{^) 
versus Q/t for the same parameters. 



m = 



- COS k'y)x"{k — k' , ri)(cos k^ — cos ky)) 
((cos kx — cos ky)'^) 

(13) 
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was calculated for U/t = 10 and Xd(^) is shown in 
Fig. ^p. These calculations are for a reduced tempera- 
ture T/Tc ~ 0.95. Once T < T^, the ti component of the 
Nambu self-energy gives (l){kA,uj) and one can calculate 
I{kA, For T/Tc = 0.95 the shift due to the magnitude 
of the gap at the antinode A(fcA) is negligible. 




FIG. 4: (a) /(/fcA,^) for the 2 x 2 DCA-NCA Hubbard cal- 
culations with different values of U/t a,t a filling (n) — 0.8 
and T/Tc ~ 0.95. (b) The density of states N{Q.) at T = Tc 
for (n) — 0.8 and different values of U /t. Here one sees the 
increase in the Mott gap that separates the upper and lower 
Hubbard bands. 

At low frequencies is seen to increase over 

the spectral range associated with the spin-fluctuation 
response seen in Xd(^)- ^ exceeds this range, /(^a, ^) 
passes through a weak maximum and dips down slightly, 
similarly to the t-J results shown in Fig. [JJa. Then how- 
ever, on a higher energy scale I{kA, ^) goes to 1. Similar 
results for I{kA, ^) for different values of U/t are shown 
in Fig. [4^. Here one sees that as U/t increases, the region 
over which I{kA,^^) remains below 1 extends to higher 
energies. For the t-J model this energy was pushed to 
infinity, but for the Hubbard model the high frequency 
contribution that takes I{kA,^) to 1 is associated with 
the upper Hubbard band as seen from the single-particle 
density of states shown in Fig. HJd. 

These numerical results show that the d-wave pair- 
ing interaction in the t-J and Hubbard models contain 



both retarded and non-retarded contributions. 241 The 



retarded contribution occurs on an energy scale which is 
small compared to the bare bandwidth 8t and the on- 
site Coulomb interaction U. For the Hubbard model, the 
"non-retarded" contribution occurs on an energy scale set 
by the Mott gap and is related to excited states involving 
the upper Hubbard band. For the t-J model, this energy 



scale is pushed to infinity and the exchange contribution 
is instantaneous. 

A simple phenomenological form for the d-wave pairing 
interaction, consistent with these observations is 

3 - 

-U^xik — k',uj,uj') — J(cos kx — cos fc„)(cos k'^ — cosk' ) . 

_ (14) 

Here co) is the dynamic spin susceptibility and U and 
J are effective coupling constants. The retarded con- 
tribution to the pairing comes from the first term, and 
the non-retarded contribution from the second, exchange, 
term. Unlike the traditional low Tc case where the non- 
retarded screened Coulomb interaction suppresses the 
gap, here the non-retarded exchange term enhances the 
d-wave gap. 

The question regarding whether there is a "pairing 
glue" is then a question of whether the dominant con- 
tribution to dAkAT^jJ = 0) comes from the integral of 
02 (^A I w)/a;.[25j From the results presented here we con- 
clude that both the t-J and Hubbard models have spin- 
fluctuation "pairing glue" . However, they also exhibit 
a smaller, non-retarded contribution. For the cuprate 
materials, the relative weight of the retarded and non- 
retarded contributions to the pairing interaction remains 
an open question. Thus the continuing experimental 
search for a pairing glue in the cuprates is important 
and will play an essential role in determining the origin 
of the high Tc pairing interaction. 
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